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ABSTRACT 



As is known, a solitary pulse in the complex cubic Ginzburg-Landau (GL) equation is 
unstable. We demonstrate that a system of two linearly coupled GL equations with gain 
and dissipation in one subsystem and pure dissipation in another produces absolutely 
stable solitons and their bound states. The problem is solved in a fully analytical form by 
means of the perturbation theory. The soliton coexists with a stable trivial state; there is 
also an unstable soliton with a smaller amplitude, which is a separatrix between the two 
stable states. This model has a direct application in nonlinear fiber optics, describing an 
Erbium-doped laser based on a dual-core fiber. 



PACS numbers: 42.81.Dp; 42.81.Qb; 52.35.Sb; 03.40.Kf 
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Localized pulses (solitons) play a central role in a number of nonlinear physical sys- 
tems, which are now a subject of a very broad interest [1,2]. Since real systems are always 
lossy, it is necessary to have an active element in the system in order to provide gain 
compensating the losses. In plasma physics and hydrodynamics, the gain is provided by 
intrinsic instabilities of the system 0. For solitons in nonlinear optical fibers (NOF's), 
an effective way to compensate for the losses is by use of Erbium-doped amplifiers [2]. 
Here, however, one encounters a fundamental problem: if the active element is uniformly 
distributed along the length of the system, it automatically renders the zero solution un- 
stable, thus lending an instability to the whole soliton. A commonly known model which 
demonstrates this property is the complex cubic Ginzburg-Landau (GL) equation, which 
in many cases, and especially in application to NOF's, may be regarded as a perturbed 
nonlinear Schrodinger (NLS) equation 0: 

iut + -Mxx + \u\'^u = ijou + i'^iu^x - i'^2\u\^u , (1) 

where u{x,t) is an envelope function (e.g., of electromagnetic waves in the NOF's), t and 
X are the spatial and temporal variables (in the NOF's their physical meaning is reversed), 
7o is the gain, while the coefficients 71 and 72 account for the dispersive and nonlinear 
losses, and all the 7's are assumed non-negative. A well-known fact is that Eq. (1) admits 
an exact solitary-pulse solution, which in the limit of a vanishing right-hand side of Eq. 
(1) goes over into the soliton of the NLS equation 0. Regarding the pulses in the model 
(1) as perturbed NLS solitons, it has been also demonstrated that they are able to form 
two-soliton and multisoliton bound states, which are stable against disturbances of the 
separation and phase differences between the solitons [Q. However, it is obvious that, due 
to 7o > 0, the trivial solution m = is unstable in this model, hence an isolated pulse 
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as a whole is unstable too (it was demonstrated numerically that locally stable pulses 
are possible in Eq. (1) with negative 70 and 72 f^; we do not consider this case here 
as the model is then globally unstable). This circumstance does not render the pulses 
meaningless objects - they may be effectively stable when the system is short enough, or 
the evolution is considered at finite times. In the general case, however, development of 
the instability leads to a dynamical chaos 0. An example is the so-called "dispersive 
chaos" experimentally observed in the binary-fiuid convection 0. 

A problem of the fundamental interest is to find sufficiently simple physical models 
which can produce totally stable pulses. A known example is the driven damped NLS 
equation |p, which has various applications, including nonlinear fiber optics [^]. However, 
in that model the pulses are not truly localized, being supported by an oscillating back- 
ground. Experimentally, absolutely stable localized pulses of traveling-wave convection 
were discovered beneath the instability onset in narrow channels filled with a binary fiuid 
heated from below |TD|. A distinctive feature of a model supporting stable pulses is bista- 



bility, as, being stable, the localized pulses must coexist with the stable trivial solution. 
The simplest possibility to provide for the bistability is to introduce a quintic GL equation 
TT| . In terms of nonlinear fiber optics, this equation models a nonlinear amplifier |jl2 



which can be built, e.g., as a combination of the linear amplifier and a saturable absorber 
with an instantaneous response. In this case, the quintic equation can be obtained by 
means of a truncated Taylor expansion of the gain/loss character isitic of the medium. 

Within the framework of the latter model, the existence of stable pulses can be demon- 
strated analytically in two opposite cases: when the quintic GL equation is close to the 
NLS equation [^, and when it is close to the real GL equation |]T^. For intermediate 



values of the parameters, existence of the stable pulses in the model was demonstrated 



numerically [jT5|. Stable localized pulses for which the quintic equation seems to be an 
appropriate phenomenological model were observed in the form of subcritical pulses in the 
traveling-wave convection in narrow channels |Ty]. Similar pulses were produced by nu- 
merical simulations of the full system of the corresponding two-dimensional hydrodynamic 
equations However, the quintic GL equation proves to be, even at the phenomeno- 



logical level, a crude model for the pulses observed in the experiment. The most essential 
feature which finds no explanation in the framework of this model is a very low velocity of 



the pulse in the laboratory reference frame |T^- As it was demonstrated by Riecke ||T8 



an essential improvement of description of the experimentally observed pulses is achieved 
within the framework of a model coupling the cubic GL equation (actually, a couple of 
such equations for counterpropagating waves) to an additional equation for a real vari- 
able, which is a "mean field" , representing in this problem a concentration field. It is well 
known that coupling of the GL equation to a real mean-field mode can drastically change 



the dynamics of the GL model [IK 



In this work, we aim to put forward a model of another type supporting absolutely 
stable solitons. Actually, the model will generate two solitons, one stable and one unstable, 
the unstable one being a separatrix between attraction domains of the stable soliton and 
of the stable trivial solution. It will also be demonstrated that the model may support 
stable bound states of solitons. 

This model describes a dual-core NOF (also known as a directional coupler ||20|), in 
which one core is active (i.e., doped with Erbium), while the other is passive (undoped). 
The couplers have attracted a lot of attention due to the possibility of applications in 
photonics In the works |22| and a coupler with one active core was proposed. 



It was demonstrated that, adding the second passive core with loose ends to the usual 



doped fiber used in the soliton-generating lasers, one can essentially improve quality of the 
generated solitons: while the soliton, being a strongly nonlinear object, remains essentially 
confined to the active core, linear noise can readily couple to the passive core, where it 
is radiated away through the loose ends. This provides effective filtering of the noise. 
Independently, it the work it was proposed to use a coupler for separating noise from 



the soliton in an existing pulse. The basic physical idea proposed in Ref. is the same 



as in Ref. ^2\: the soliton keeps itself in the core in which it is propagating, while the 



noise easily couples to the second core. It was assumed in Ref. pj| that the second core 



was lossy, so that the noise tunneling to that core would be killed there by the losses. An 



important result obtained both analytically and numerically in Ref. is that the best 
efficiency of the filtering would be attained not at very strong losses in the second core, 
but in the case when the loss coefficient is close to the coupling constant between the two 
cores. 



The numerical simulations reported in Refs. - EM demonstrate that the soliton 



seems very stable in these models (moreover, in Ref. it was demonstrated that the 
model produced a robust soliton even in the case when dispersion in the active core was 
normal, i.e., the solitons could not exist in the usual NLS equation). These observations 
suggest to analyze existence and stability of solitary pulses in the two-component models. 
A remarkable fact is that, as it will be shown below, this problem can be consistently 
solved in a fully analytical form. 

We adopt the following model, which is a combination of the more special ones from 
Refs. [^] and ||2^ (unlike Eq. (1), here we will use the standard "optical" notation for 
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the spatial and temporal variables z and r): 

+ l^TT + \U^U — i'JoU — i'^lUrr + Kt> = 0, (2) 
iVz + ]^VrT +\v\^V + iTqV + KU = (3) 

where the variables u and v are envelopes of the electromagnetic waves in the active 
and passive cores, the coefficients 70 and 71 are the same as in Eq. (1), the nonlinear 
dissipation is neglected (it can be easily restored), k is the coupling constant, and Fq is 
the loss coefficient in the passive core. An additional loss coefficient Fi in the passive core, 
similar to 71 in the active one, may be added to the model. However, the additional lossy 
term would render the following analysis more cumbersome without making any essential 
change. Contrary to this, it will be demonstrated that keeping the term ~ 71 in Eq. (2) 
is necessary. 

As concerns physical applications of this model, it was already mentioned that it was 
closely related to the ones describing a soliton-generating laser |Q (but in that case, 
Fq = 0) and a time-domain fiber ffiter 0] (however, 70 = 71 = in the latter model). 
Very recently, it was demonstrated that the model based on Eqs. (2) and (3) in their 
full form finds another practically important application: it is a basis for a new type of a 
nonlinear optical amplifier 



First of all, we will consider stability of the trivial solution u = v = 0. Inserting 
into linearized Eqs. (2) and (3) {u^v) ~ ex\i{az — iuit), one can immediately obtain the 
relation between the instability growth rate a and the perturbation frequency u: 

+ (Fo - 7o + 7i^^^) a + k'^ - Fo7o + 'JiTquj'^ = 0, a = a+ ^iuj"^. (4) 

The stability condition for the trivial solution implies that Re a < at all real u, which is 
equivalent to demanding that the free term and the coefficient in front of the linear term 



in Eq. (4) are always > 0. Eventually, this amounts to 

7o < To < K^ho. (5) 

Obviously, a necessary condition following from Eq. (5) is 70 < i.e., the coupling must 
be stronger than the gain. Notice that the coefficient 71 does not appear in Eq. (5). 

Now, let us consider evolution of the soliton in the system of Eqs. (2) and (3). We 
will treat this problem perturbatively, assuming that the coupling, gain, and losses are 
all small perturbations to the NLS equation, although different perturbations may have 
different orders of smallness. Actually, we will assume that the gain and losses in the 
active core are essentially weaker than the coupling between the cores, while the losses in 
the passive one may be comparable to the coupling. Anyway, the perturbative treatment 
of all these terms is quite reasonable in application to the NOF's. 

The soliton has two nontrivial parameters, viz., the amplitude and velocity. The 
simplest way to derive perturbation-induced evolution equations for these parameters is 
to use the so-called balance equations for the energy and momentum of the soliton [1]. 
In our model, the evolution equation for the velocity will be exactly the same as in the 
model (1), which describes the process of "braking" of the soliton by the linear frequency- 
dependent losses [1] , therefore we will not consider this equation, and will simply set the 
velocity equal to zero. Then, in the zeroth-order approximation, the soliton resides only 
in the first core and has the form 

= 77sech(7^T) e*'*^^^\ (6) 

where ^ = ^Tf, and r} is the soliton's amplitude. In the next approximation, one seeks for 
a component ("shadow") of the soliton in the second core [1]. Obviously, it has the form 

u{z,t) = 1^(t) exp(i0(^)), where the real function V{t) is determined by the equation 
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following from Eqs. (3) and (6): 



rfV = —2K,r]sech{r]T) (7) 



(at this step, we neglect the lossy term in Eq. (3) which enters only at the next order). 
Finally, the slow evolution of the soliton's amplitude under the action of the perturbations 
is determined by the above-mentioned balance equation for the quantity which plays the 
role of energy in nonlinear optics: 

/+00 , , 
(\u{r)\^ + \v{r)\'')dT. (8) 

An exact balance equation following from Eqs. (2) and (3) is 

^ = 270/ \u{r)\'dr-2^J \ur(r)\'dr-2ro \v(r)\'dr. (9) 

dz J —00 J— 00 J —00 

Now, one should find V{t) from Eq. (7) and insert it into the last term on the right- 
hand side of Eq. (9), while a contribution from v to the left-hand side (see Eq. (8)) may 
be neglected in the lowest nontrivial approximation. Eq. (7) can be solved by means of 
the Fourier transformation, which leads to the following integral representation for V{t): 

sech {tit') e"''!^-^ ^dr'. (10) 

-00 

Inserting this expression into the last term of Eq. (9), one can explicitly calculate all the 
integrals, which eventually yields the evolution equation sought for: 

^ = 27077 - l^iv' - Ck'ToV-', (11) 

where C = |7r^ + C(3) ^ 2.845. The formal singularity of the last term in Eq. (11) at 
77 ^ is fictitious, as this expression is irrelevant at very small 77, see below. 

It is straightforward to see that Eq. (11) gives rise to two physical {rj'^ > 0) fixed 
points, provided that 

7o' > {3C/8)k'To^I, (12) 
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and to no fixed points in the opposite case. Thus, (12) is the necessary and sufficient 
condition for existence of the sohtons in the considered modeL It is, of course, important 
to check if this condition is compatible with the other fundamental condition, (5), which 
is necessary for stabihty of solitons in the model. Because Eq. (5) does not involve the 
parameter 71, one can secure the compatibility simply choosing 71 to be small enough. 

Next, it is easy to check that, once the condition (12) is met, the solution with larger 
rf is stable, and the one with smaller rf is unstable. It is very plausible that the sohton 
corresponding to the larger root is a completely stable solution in the full model, while 
the smaller root corresponds to an unstable soliton which plays the role of a separatrix 
between the stable sohton and the stable trivial solution. 

Usually, existence of solitons is related to modulational instability of continuous wave 
(cw) solutions [2]. The cw modulational (in)stabihty in the present model will be consid- 
ered in detail elsewhere. However, in the regime which is akin to the case considered in 
this work, i.e., when the coefficients of the gain and loss are small, the coupling constant 
is small too, and the field in one core {v) is therefore much weaker than in the other {u), 
the modulational instability in the present model is, evidently, close to that in the single 
NLS equation. Thus, as well as in the usual NLS equation, one may qualitatively regard 
the solitons in the present model as pulses produced by the modulational instability of 
the cw. 

It is now relevant to discuss conditions guaranteeing application of the perturbation 
theory to this problem. The "primary" conditions are 71 <^ 1, \n\ -C rf., and Tq <^ rf. 
Since, in a typical case, rf ~ 70/71, the final set of the applicability conditions takes the 
form 

7i -C 1, 7iro -C 7o, 7i|«;| < 7o- (13) 
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Obviously, these conditions are compatible with the underlying inequalities (5) and (12). 

When formulating the model, we have omitted the dispersive loss term iTiUrr in Eq- 
(3). If kept, it will generate an additional term ~ —k'^TiI]^^ in Eq. (11). It is easy 
to check that the latter term will not produce any qualitative difference in properties of 
the fixed points. On the other hand, if the term ~ 71 is omitted in Eq. (2), the result 
will be disastrous: there will remain a single unstable fixed point. Actually, the model is 
globally unstable in that case. At last, if the nonlinear lossy term ~ 72 (the two-photon 
absorption, in terms of the nonlinear optics) is added to Eq. (2) (cf. Eq. (1)), it will 
merely renormalize the coefficient 71 in the final results displayed above. 

Now, we will briefly consider bound states (BS's) of the solitons in this model, fol- 
lowing the lines of Ref. The BS's may exist due to the fact that the small linear 
terms accounting for the gain and dissipation render solitons' tails oscillatory, which, in 
turn, gives rise to local minima in the effective potential of the soliton-soliton interaction, 
produced by overlapping of the "head" of each soliton with the tail of the other one. How- 
ever, in the framework of the perturbation theory, the BS's are fragile, although stable: 
the distance between the solitons in the BS is large, and, accordingly, the corresponding 
binding energy is exponentially small. Nevertheless, existence and stability of the BS's 
predicted by the perturbation theory was confirmed by direct numerical simulations |2^ 
of the driven damped NLS equation (see also 0); recently, this prediction was also 



confirmed, with a fairly good accuracy, for the cubic GL equation (1) In the present 
model, the BS's can be rendered more robust by increasing the dissipative constant Fq in 
the passive core. Therefore, we will consider the case Fq ^ 70, 71?]^, which is compatible 
with all the above conditions necessary for existence of the stable solitons. In this case, 
we again consider linearized Eqs. (2) and (3); however, instead of the plane- wave solution 

11 



leading to the dispersion equation (4), we are now interested in an exponentially decaying 
solution describing the soliton's tail: 

{u,v) ~ exp(— ?7|r| + ix|r| + zgz) , (14) 

where we consider rj as a given soliton's inverse size (cf. Eq. (6)), while x ^i-nd q must be 
found. The linearized equations immediately yield (with regard to the assumed dominance 
of To): 

X = ro/2r/, (g - + Ix')' = n'- \tI. (15) 

According to Ref . , the minimum separation T between solitons in the BS is determined 
by the coefficent x Eq. (14): T = ti/2x = tttj/Tq, where we have made use of Eq. (14). 
On the other hand, the second relation in Eq. (15) imposes a fundamental limitation 
Fq < otherwise the soliton simply does not exist. This, in turn, leads to a limitation 
on the minimum separation between the bound solitons: T > TTri/2K,. 

In conclusion, we have demonstrated that a simple analytically tractable model, based 
on the linearly coupled cubic GL equations, admits a fully stable soliton coexisting with 
the stable trivial state. The model finds a direct physical realization in terms of the fiber 
laser and suggests a way to stabilize soliton generation in the laser. 

This research was partially supported by the National Science Foundation through the 
Center for Ultrafast Optical Science at the University of Michigan, Ann Arbor. 
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